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The neutrino flux close to a supernova core contributes substantiaUy to neutrino refraction so that 
flavor oscillations become a nonlinear phenomenon. One unexpected consequence is efficient flavor 
transformation for anti-neutrinos in a region where only neutrinos encounter an MSW resonance 
or vice versa. Contrary to previous studies we find that in the neutrino-driven wind the electron 
fraction Ye always stays below 0.5, corresponding to a neutron-rich environment as required by 
r-process nucleosynthesis. The relevant range of masses and mixing angles includes the region 
indicated by LSND, but not the atmospheric or solar oscillation parameters. 
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Introduction. — The evidence for flavor oscillations of 
solar and atmospheric neutrinos and in the LSND exper- 
iment implies mass differences so small that refractive 
effects influence or even dominate neutrino oscillations 
in many situations of practical interest. However, there 
are only two examples where neutrinos themselves as a 
medium modify the oscillations. One is the early uni- 
verse , the other core-collapse supernovae (SNe) |2| . 

In a seminal paper Pantaleone showed that neutri- 
nos as a background medium differ markedly from other 
fermions. A given background neutrino may be a co- 
herent superposition of flavor states, implying an "off- 
diagonal refractive index" in flavor space. The oscil- 
lations of the entire ensemble thus become a nonlinear 
phenomenon with unexpected consequences. When the 
neutrinos themselves dominate as a background medium, 
the oscillations become "synchronized," i.e. all modes 
oscillate collectively with the same frequency, a behav- 
ior first discovered by Samuel With our collabora- 
tors we recently found a simple physical interpretation 
of this perplexing phenomenon in terms of the dipole- 
dipole coupling of a collection of magnetic dipolcs which 
spin-precess in an external magnetic field Q . 

The first environment where v-v refraction plays a 
crucial role is the epoch of the early universe that 
precedes big-bang nucleosynthesis. If initially large 
flavor-dependent i^-i?-asymmetries exist, they may be 
equilibrated by oscillations and collisions before weak- 
interaction freeze-out so that the primordial helium 
abundance implies stringent limits on the overall cos- 
mic neutrino density 0, |8[ ^. Depending on initial 
conditions the modiflcation of the flavor relaxation pro- 
cess caused by the synchronization effect is only mild, 
or it may even prevent equilibrium entirely because the 
synchronized oscillation frequency can become arbitrarily 
small. The interplay of simultaneous v and i> oscillations 
is a crucial and non-trivial ingredient in the evolution of 
this system. 

The second system where background neutrinos may 
be important is the rarefled region just outside the 



nascent neutron star a few seconds after SN core bounce. 
In the innermost regions of the SN, the neutrino flux 
is so large that the weak-interaction potential created 
by the neutrinos is comparable to that of the ordinary 
medium. The neutrino spectra and fluxes differ between 
the flavors and between neutrinos and anti-neutrinos of 
a given flavor. Swapping the fluxes of different flavors 
would crucially modify the production of heavy elements 
via r-process nucleosynthesis if th is p henomenon takes 
place in the SN hot bubble region |l^ . The possible 
importance of i/-i/-refraction in this context was quickly 
recognized jlj, [ij, 0. The main consequence implied 
by these approximate treatments was a small shift of the 
oscillation parameters where a significant spectral swap- 
ping by resonant oscillations takes place. 

Alerted by the subtleties encountered in our study of 
early-universe oscillations Q we revisit the v-v-eSect in 
the SN hot bubble region. We find that previous au- 
thors indeed underestimated the complications that arise, 
in particular, when neutrinos and anti-neutrinos oscil- 
late simultaneously and cause refractive effects for each 
other. We find, for example, that in a region of param- 
eters where neutrinos encounter an MSW resonance, the 
anti-neutrinos are "dragged along" and also show large 
flavor transformations. The final picture of the interplay 
between neutrino oscillations and r-process nucleosynthe- 
sis is very different than previously imagined. 

Two-Flavor System. — To be specific we study the Ve-Vii 
system with oscillation parameters tan^ 9 and Am^ = 
m\ — m\ > 0. In the absence of neutrino background ef- 
fects, neutrinos (anti-neutrinos) encounter an MSW res- 
onance for tan^ 6 <1 (tan^ 9 > 1). The evolution of the 
neutrino system is described by the 2x2 density matrices 

Pv{t)^{^:' '';']^\[p,{v,t)+cT-v^{t)] , (1) 

and analogously pp for anti-neutrinos. Here, Ci are the 
Pauli matrices while Pp(i) and Pp(i) are the usual po- 
larization vectors for v and v modes with momentum 
p, respectively. The diagonal elements paa(^^t) are the 
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occupation numbers of flavor a with momentum p. 

In the region of interest neutrinos stream freely so that 
we may ignore collisions. Therefore, the radial evolution 
equation is the usual precession formula, augmented by 
the v-v refractive term 



dqCpq(Pq-Pq 





(2) 



Here B = (sin 26*, 0, - cos 26*) is a "magnetic field," B the 
vacuum mixing angle, and z a unit vector in the z-direc- 
tion in flavor space. Further, = Y^Nb is the electron 
density with the electron fraction and Ng the baryon 
density. Finally, Cpq = 1 — p • q, implying that coUinear 
neutrinos do not cause a mutual refraction effect. 

As a matter density profile for the hot bubble region we 
use the one shown in Ref. jl^ which roughly falls off as 
r~^. As a boundary condition we assume equal luminosi- 
ties Li, for all flavors of order Lq = 10^^ erg s~^. The 
spectra are taken to be Fermi-Dirac distributions with 
mean energies {E^J = 11 MeV, {Ej^J = 16 MeV and 
{^v^.i>^) — 25 MeV, respectively. These choices may not 
be entirely realistic [p^ , but for consistency with previ- 
ous work we stick with these traditional assumptions. 

In the absence of oscillations and for radially moving 
neutrinos the diagonal elements of the density matrix at 
radius r and for neutrino momentum p are 
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P 





47rr2 77r4 T^^ exp(p/T^„) + 1 



(3) 



where T^^ — (E^^) /3.151. For the iz-i' refractive effect 
the angular divergence of the neutrinos is crucial. As in 
previous works Q we use a flux-averaged value, i.e. 
in Eq. (0) we substitute 



dqCpq(Pq - Pq) X Pp ^ F (r) (P - P) X Pp . (4) 



Here, P and P are the total polarization vectors and 
F{r) = i[l — (1 — i?^/r^)^/^] is a geometrical factor with 
Ri, the neutrino-sphere radius (see for a more de- 
tailed discussion of the geometrical dependence). Both 
F{r) and the luminosity fall off as so that the v-iy 
refractive term scales as r^* at large r. In the neutrino- 
driven wind phase the medium density typically falls off 
as so that at large distances the ordinary medium 
dominates. However, at distances of 15-30 km the neu- 
trinos may dominate. 

R-Process Nucleosynthesis. — A key necessary condi- 
tion for this process to occur in the SN hot bubble 
region is that the environment must be neutron-rich. 
The neutron-to-proton ratio is fixed by the /3 processes 



+ n ^ p + e~ and Ve + p n + while charge neu- 
trality requires n/p = l/Y^ ~ 1 JTo| . Therefore, a min- 
imal requirement is Yg < 0.5, but a successful r- process 
may require Yg < 0.45. Near weak-interaction freeze-out 
(WFO), at a radius 30-35 km, only the direct (3 processes 
are important and the electron fraction is 



(5) 



where e = {E^ )/{E^^), and I the analogue for Ve- We 
take the Ve and Ve cross sections on nucleons to be equal, 
see however 
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In the absence of neutrino oscilla- 
tions and with our choice of neutrino flux parameters one 
finds Ye — 0.41, allowing for a successful r-process. 

Spectral Swapping by Oscillations. — If neutrino oscilla- 
tions occur within the WFO radius, the effective i^^ and 
Df, flux spectra change and modify Y^. As a first example 
we use Am^ = 10 eV^ and tan^ 6 — 10^^ which yield the 
Ye profile shown in Fig. |l|. For the curve marked 0, neu- 
trino background effects were ignored, the other curves 
are for the indicated values of L^. 

The oscillations can be calculated analytically in the 
limit 3> Lq where the neutrino background strongly 
dominates. We define I = P — P, integrate Eq. (||) over 
the neutrino spectra to get the evolution equations for P 
and P, and subtract them to obtain 



9,1 = dp 



~2i 



B X [Pp + Pp] + \/2 GpiVe z X I . (6) 



The neutrino background term is proportional to Ix I and 
thus vanishes. However, the individual modes Pp and Pp 
precess fast around I as in Ref. . The evolution of I is 
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FIG. 1: Ye as a function of radius for the indicated choice 
of oscillation parameters. The labels indicate Li, in units of 
Z/o = 10^^ erg s~^; — Q implies the absence of neutrino 
background effects. 
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No V background = 0.1 Lo L, = L.-, 




bg(Lan^'i^) log(ton^i>) log(tan^'i>) 



FIG. 2: Spectral swapping as a function of Am^ and X.sxi' Q for — Q (no neutrino background effects), 0.1, and ILq. The 
solid contours indicate the LSND+KARMEN2 allowed region while the region inside the dashed contour is excluded by 
the Bugey experiment [pO[. Upper Panels: Ye at the WFO radius of ~ 30 km. Lower Panels: L{ue)/L„ at r = 50 km. 



a slow precession with a certain synchronized frequency 
^synch- Wc express Wsynch by the neutrino momentum 
Psynch that would prcccss with LUgynch in the absence of 
a neutrino background. Note that the synchronization 
of both neutrino and anti-neutrino modes occurs despite 
the presence of a CP asymmetric background 

To find Psynch we use that in the present limit all 
Pp and Pp are essentially aligned with I so that their 
projections along I are conserved. All modes start in 



the z-direction so that altogether Pp 
Pp ~ Pz{p) I- We can then rewrite Eq. 

drl^ B + V2GFiVeZ 



P,(p)I and 

as 



where 



Psynch 



2psynch 
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X I. 



T: 



(7) 



(8) 



We have used that initially T^^ — Tjj^^ . For our assumed 
spectra we find Psynch — 2.2 MeV, much smaller than typ- 
ical energies of the neutrino spectra. As a consequence. 



neutrino oscillations are effective at smaller radii and for 
smaller Am? than without a neutrino background, an 
effect already observed in Refs. [^3[ 

The results of Fig. |^ are now easily explained in two 
limiting cases. Without neutrino background all neutrino 
modes experience an independent adiabatic MSW tran- 
sition starting at low energies (Ye decreases) until the 
entire z/e spectrum is swapped with that of z/^ , leading to 
the asymptotic value ~ {I + {E^^) / {El^))^'^ ~ 0.61. 

The other limiting case (L^ ^ Lo) corresponds to a 
synchronized MSW transition of the entire neutrino and 
anti-neutrino ensemble, where all modes follow an adi- 
abatic transition at the same radius where a neutrino 
with momentum Psynch would do an MSW transition in 
the absence of background neutrinos. Y^ takes on the 
value 0.5 because both and are swapped with 
and v^, respectively, and thus take on identical spectra. 

For the intermediate cases there is some degree of syn- 
chronization, but it is gradually lost at larger radii with 
the dilution of the neutrino flux. Still for the nominal 
neutrino luminosity = Lq the evolution is quite dif- 
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ferent from the no-background case. 

For our assumed flux spectra we have systematically 
calculated the effect of spectral swapping as a function of 
Am^ and tan^ 9. In Fig. || (for black-and-white printing 
see Fig. ^ we show our results for the assumed luminosi- 
ties = (no neutrino background effects), L^, ~ 0.1, 
and 1L{). We indicate the region of mixing parame- 
ters which is compatible with the experimental results 
of LSND and KARMEN2 from a joint analysis |9| and 
the region excluded by Bugey |20[ . 

In the upper panels we show at the WFO radius 
~ 30 km. In the absence of neutrino background effects 
{Ln = 0) our results agree with those from the previous 
literature . For instance large Am^ and small tan^ 9 
cause Yf, > 0.5, violating the minimal requirement for r- 
process nucleosynthesis. However, such regions gradually 
disappear when the neutrino background is enhanced, 
i.e. neutrino background effects prevent Ye from exceed- 
ing 0.5, in stark contrast to the previous literature. Like- 
wise, for an inverted mass situation (tan^ 9 > 1) spectral 
swapping effects are quite significant even though neutri- 
nos do not encounter an MSW resonance. 

Previous studies of the neutrino background effect used 
various approximations ||l2, 13, lj|. Rcfs. |lj] consid- 
ered the full set of equations in an approximate way, but 
did not include the evolution of anti-neutrinos. We in- 
terpret the fundamental difference between the previous 
literature and our results as being caused by the simul- 
taneous oscillations of neutrinos and anti-neutrinos. 

Another intriguing aspect of our results is the behav- 
ior of Livf.) which we show in the lower panels of Fig. ^ 
for r = 50 km. It is approximately this region where 
the protons and neutrons of the material which would 
eventually undergo heavy-element synthesis form alpha 
particles |2|, ^ Therefore, an equal number of 

protons and neutrons is locked into alphas, the excess 
of either of them remaining free. But an excess of neu- 
trons can be erased by capture if L(Pf,) is large enough 
("a-effect"). One speculative way of reducing L{i'e) m- 



electron fraction Yg was never enhanced above 0.5 when 
neutrino background effects were included, thus fulfilling 
the minimal condition for r-process nucleosynthesis. 

If the LSND signature is not due to neutrino conver- 
sions and the active-active oscillation parameters perma- 
nently settle in the regions indicated by solar and at- 
mospheric neutrino conversions, then it is unlikely that 
neutrino oscillations infiucncc r-process nucleosynthesis, 
always assuming the SN hot bubble region is the cor- 
rect site. There remain interesting oscillation effects at 
larger radii where the matter density is smaller , but 
we do not expect ^'-^-refraction to play a major role at 
these distances. This caveat notwithstanding we find the 
nonlinear effects of i/-z^-refraction in the SN hot bubble 
region a fascinating topic worth investigating. 
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